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We discuss recent developments in multiloop calculations aiming eventually in computing the total cross section 
for e + e~ annihilation into hadrons a^ad in order O(a^). 
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1. Introduction 

Grand Unification, the merging of the three 
gauge theories based on f7(l), SU(2) and SU(3) 
groups with their three independent coupling g\, 
<?2 and gs into the unified framework of SU(5) 
or, more probable, 50(10) gauge theory is one 
of the most attractive possibilities for physics be- 
yond the Standard Model. As is well known, this 
can only be considered as attractive and natural 
option consistent with the present values of the 
coupling constants if supersymmetry is realized in 
Nature. Indeed, most of the models predict that 
at least some supersymmetric partners of quarks, 
leptons, gauge and Higgs bosons will be discov- 
ered and studied in detail at the next generation 
of colliders, LHS or ILC. 

Once the masses of sufficiently many SUSY- 
partners are measured precisely, the detailed en- 
ergy dependence of the coupling constants will be 
fixed and the test for the unification can proceed 
with significantly improved precision. In fact, it 
is conceivable that some hints about symmetry 
breaking mechanism responsible for the transi- 
tion from, say, 50(10) to SU(1) x SU(2) x 5f7(3) 
could be inferred from such an analysis, in partic- 
ular, if combined with a similar study for Yukawa 
unification. 

However, for such a program to be pursued 
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successfully, the initial values entering into the 
evolution equations — gauge couplings and masses 
as measured at low energies — have to be known 
with sufficient precision. Of particular interest is 
the gain in resolution power from a more precise 
determination of the strong coupling. A signifi- 
cant improvement is anticipated for the GIGA- 
Z option of the linear collider. Indeed, it has 
been argued in P that a sample of 10 9 Z decays 
might well lead to an experimental uncertainty of 
Sa s — 0.0008, and even a value of Sa s — 0.0005 
has been quoted [2]. 

The determination of a s from Z is either based 
on the measurement of the leptonic branching ra- 
tio 

Ri = T h /T e = 20.767 ±0.025 

which amounts to a simple counting of leptonic 
and hadronic final states, or it exploits the lep- 
tonic peak cross-section 



r 2 

i tot 



= 2.003 ± 0.0027 nb. 



In the Standard Model fit 
a s = 0.1183 ±0.0030 

all this information is combined py. In all these 
cases the a s dependence enters through the for- 
mula 
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where the second factor, conventually denoted 
as l/3R(s), gives the perturbative series valid 
for massless quarks. The "corrections" in the 
above formula originate from quark mass terms 

2 

~ j-rt, from singlet terms which start in order 

a 2 s , which present in the axial amplitude only and 
which originate from the imbalance between top 
and bottom quark loops and, finally, mixed, non- 
factorizeable QCD x electroweak corrections of 
order a s a woak . 

All these terms are discussed in 0], and, indeed, 
the "corrections" terms in eq. are typically 
well under control. The dominant theory error 
one encounters in the extraction of a s originates 
from estimates of the not yet available terms of 
order a A s . To estimate this uncertainty, different 
strategies have been advocated and we shall only 
list a small subset of the possible choices. 

The most conservative estimate of the trunca- 
tion error for an asymptotic series is identical to 
the last calculated term, presently of order 
and leads to an uncertainty of Sa s — 0.002 cor- 
responding to 5a s /a s = 1.8%. Alternatively one 
may vary the renormalization scale which appears 
as argument of the strong coupling /j, in a plau- 
sible range around y/s , say from 1/3-s/s to 
Using the correspondingly modified perturbative 
series to estimate a s (//) (and evolving a s (p) back 
to a s (Mz)) one obtains a quite asymmetric varia- 
tion Sa s = lo'ooo2 °f the same order as estimated 
above. 

Various prescriptions have been used to esti- 
mate not yet calculated higher order terms. From 
the Principe of Minimal Sensitivity (PMS) 
or Fastest Apparent Convergence (FAC) [Hj the 
forth order term is predicted as — 97(a s /-7r) 4 [7|. 
This then would lead to a shift in a s of 0.0006, 
corresponding to 8a s /a s = 0.5%. 

From this considerations it is evident that for a 
reliable determinations of a s from GIGA-Z, fully 
exploiting the experimental precision, the evalu- 
ation of the a 4 , term is a must. This arguments 
get even stronger in those cases, where low energy 
measurements are involved. The most extreme 
case in this direction is the determination of a s 
from the semileptonic versus leptonic decay rate 
of the tau lepton Considering the enormous 



events rates at low energy electron-positron stor- 
age rings, like B- or charm- or Tau- factories, R 
measurements at 10.5 GeV or at 3.7 GeV, just be- 
low the threshold for open bottom or charm pro- 
duction, could also lead to fairly precise measure- 
ments of a s . Indeed, to achieve the benchmark 
precision Indeed, to achieve the benchmark preci- 
sion 8a s {Mz) — 0.003 would require Sa s = .007 
at 10.5 GeV and 6a s = .013 at 3.7 GeV Evidently 
only systematic uncertainties are the limiting fac- 
tors for an ultra precise determination of a s . 

2. The long march towards R in a 4 

Let us now briefly describe the strategy to eval- 
uate the a 4 term and the status of this calculation 
and a variety of physically and mathematically 
relevant results that are available at present. 

Consider the correlator of two currents (all 
Lorentz and flavour indices are suppressed and 
quark masses are set to zero) j = ipTip and 
ji = tpT^tp 

rP'V) = * / dxe^{Q\T\j(x)3\ti)W) 

which is related to the corresponding absorptive 
part R(s) through 

R jj (s) KSff' J '(s -iS). 

Rcnormalized and bare correlators are related 
through 

II- z ! - ■ ■/.• n;;; 7 -.,»^i 

and the independence of II s on the renormaliza- 
tion scale is reflected in the renormalization group 

2 

equation of the form [L = In 3^2) 

(JL _ 2 y( Qa ) +p(a s )-^j U" = j"(a s ). 

(2) 

Here the anomalous dimension 7 J (a s ) is related 
to the corresponding renormalization constant as 
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while the QCD beta-function 

Equation J2J clearly demonstrates that in order 
to find the q 2 dependent part of II at (N+l)- 
loop (corresponding to af order) one needs the 
(N+l)-loop anomalous dimension and the N- 
loop approximation of II, (i.e. of order a^ 1 ) 
including its constant part. Note that the oper- 
ation ct s -^- raises the power of a s by one unit 
as (3(a s ) starts from a 2 . The problem of finding 
of a (N+l)-loop anomalous dimension, however, 
can be reduced in a systematic automatized way 
to the evaluation a proper combination of N-loop 
massless propagator integrals |9I10| . 

The situation is significantly simplified when- 
ever ~f n happens to vanish. One particular and 
physically relevant case is the m 2 /s term in the 
small mass expansion of the absorptive part of 
the vector correlator, where the corresponding 
four-loop 0{aP s ) contribution to II yy is indeed 
sufficient to obtain the term of order 0{a i s m 2 q / s) 



2.1. Reduction of massless propagators 

Thus, to obtain the massless i?-ratio in order 
a.g, a large number of four-loop massless propa- 
gators, including their finite parts, must be eval- 
uated. The complications can be best judged by 
contrasting the cases of the three- and four-loop 
calculations: 11 (about 150) topologies with and 
3 (11) topologies without insertions are involved 
in the three-(four-)loop case. The reduction has 
to be performed to 6 (28) master integrals. Most 
importantly, however, a set of recursion relations 
based on integration by parts identities is avail- 
able for the 3-loop case. These recursion relations 
have been constructed manually ^Sj and imple- 
mented in the program MINCER 

A straightforward implementation of this con- 
cept to the four-loop case seems to be difficult at 
present. An alternative concept has been advo- 
cated in 15 which allows to perform the reduc- 
tion of an arbitrary massless propagator to a sum 
of master integrals "mechanically" in the limit of 
large space-time dimension d. Consider the re- 
duction of an amplitude f(d) which depends on 



the topology, the power of the various propaga- 
tors and the space-time dimension d to master 
integrals: 

f(d) = J2 C a (d)*M a (d), 
0: = masters 

where M a , with a = 1 — 28 stand for master 
integrals. The coefficients C^ a (d) are known to 
be rational functions of d; their determination is 
the central problem to be solved. In the limit of 
large d we have 

h 

The terms in the 1/d expansion can be expressed 
(with the use of the new representation for Feyn- 
man integrals developed in |lfill7| ) through sim- 
ple Gaussian integrals — a task of purely algebraic 
nature. Obviously, given sufficiently many coef- 
ficients, the functions C a can be reconstructed. 
However, the degrees m and n are strongly de- 
pendent on the power of the propagators in the 
integral under consideration with drastic effects 
on the effort involved in the evaluation. In the 
process of reduction a 4-loop diagram of a typical 
topology to Gaussian integrals one should handle 
a polynomial of 9 variables of degree k consist- 
ing of terms. For k=24 this corresponds 
to approximately 4 • 10 7 terms or 4 GB storage 
space, while k=40 leads to 2 • 10 9 terms, corre- 
sponding to 200 GB. Weeks, partly even months 
of runtime and hundreds of GB disk space are re- 
quired for the evaluation. On the other hand, the 
purely algebraic manipulations required to per- 
form the Gaussian integrals are ideally suited for 
treatment within the algebraic manipulation pro- 
gram FORM [TBI and its parallel version PAR- 
FORM |19I20| . which both are specially tuned for 
dealing with a huge volume of algebraic manip- 
ulations. The method of large d expansion has 
been successfully applied for solving a number of 
problems which will be briefly discussed in the 
next section. 

2.2. Important dots 

The reduction of the (N + 1) loop anoma- 
lous dimension to the calculation of A-loop mass- 
less propagators discussed above works "natively" 
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only for logarithmically divergent integrals while 
H n is in general quadratically divergent. The only 
known way to apply the reduction in such cases 
is to use (double!) differentiation w.r.t. the ex- 
ternal momentum q to decrease the dimension. 
This lead to " dots" = squared propagators which 
immensely complicate all calculations. 

An important and non-trivial simplification ex- 
ists for the scalar (SS) correlator due to the well- 
known Ward identity: 

%9^JiM) = K T m i ) 2 nf/ P (5) + OK). 

Basically this means that the 0(m 2 ) part of the 
longitudinal part of VV correlator is identical to 
the massless SS one. This allows to compute 
the 0(rriq) part of the VV correlator instead of 
the massless SS which resulted in diagrams with 
one squared propagator less and saves a lot 
of work! Unfortunately, no such simplification is 
known for the case of the massless vector correla- 
tor. 

3. Results 

Although the R ratio to order af for the vec- 
tor correlator is not yet available a number of in- 
termediate results of phenomenological relevance 
have been obtained recently. 

• All 28 master integrals have been evaluated an- 
alytically ini- 

• The terms of order a A nj and a\nj for the R 
ratio discussed in section 1 were obtained in |22l 
I23p (a s = a s /-K, dots stand for not yet computed 
term of order n/ and nj) 

R(s) = 3 jl + a s + a 2 s (1.986 - 0.1153 n f ) 

+ a 3 (-6.6369 - 1.2001 n f - 0.00518 nj) 
+ a 4 (0.02152 nj - 0.7974 n} + ...)\. 

• The full dependence of the r-lepton decay rate 
on the strange quark mass m s up to order a 3 has 
been obtained in |24| . 

• The quadratic term in the small quark mass 
expansion of the i?-ratio 

i?( S ) = 3 {^ + ^} + .... 



has been obtained in ^2]- For the number of ac- 
tive quark flavours nf = 4 the result reads 

r^/12 = a s + 9.09722a 2 

+ 52.913a 3 , + 128.499 a% (3) 

Let us mention that the a 4 term is predicted to 
177 and 193 by the estimates based on PMS and 
FAC respectively. 

• Most recently the five loop anomalous dimen- 
sion of the scalar correlator was obtained 25 . In 
combination with the finite part of the four-loop 
scalar correlator this corresponds to the R-ratio 
of the scalar correlator (denoted with an extra 
tilde below) and has important applications for 
QCD sum rules and for the Higgs decay rate to 
b-quarks |25l26j . Written for brevity in numerical 
form only, the result reads: 

R = l + 5.6667a s + [35.94- 1.359 n f ] a 2 s 

+a 3 [164.1 - 25.77 n f + 0.259 nj] + (4) 

a 4 s [39.34 - 220.9 n f + 9.685 nj - 0.02046 nj] . 

In Table QJ the result for r± - the coefficient in 
front of a* term in (0} - is compared with pre- 
dictions obtained in works |27I28| . Note that the 
Principle of Minimal Sensitivity or the Principle 
of Fastest Apparent Convergence used in pro- 
duce identical result at order a A s . The two predic- 
tions of FAC /PMS for f 4 correspond to either the 
consequence of the prediction for the correspond- 
ing Euclidean quantity (second line) or to the di- 
rect application of FAC /PMS to estimate f 4 (the 
third line). As a consequence of the large can- 
cellations in r4 the second prediction looks much 
better than the first, despite the fact that the 
estimation of the corresponding Euclidian coeffi- 
cient is quite close (within 10%) to the exact re- 
sult (for more details see |2SI)- The Asymptotic 
Pade-Approximant Method (APAM) estimation 
of r4 constructed in |2Hj fails to reproduce even 
the sign of the exact result. Finally, predictions 
of the prescription proposed by Brodsky, Lepage 
and Mackenzie (BLM) 29^ for the n/ dependent 
terms of order a A s have been communicated to the 
authors 3 : af(-260n/ + 13 7^-0.046^) and are 
also in reasonable agreement with the exact result 
of eq. (0J. 

3 M. Binger and S. Brodsky, private communication. 
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Table 1 

Comparison of the results for with earlier esti- 
mates based on PMS, FAC and APAM. 
nf 3 4 5 



fi (exact) 

h (EZI, PMS, FAC) 
h (E3, PMS, FAC) 
r 4 (E3, APAM) 



-536.8 -690.7 -825.7 

-945 -1099 -1237 

-528 -749 -949 

252 147 64 



4. Conclusion 

The calculation of the i?-ratio of the scalar cor- 
relator at order a A s has clearly demonstrated the 
enormous complexity inherent into any calcula- 
tion of such multiloop level in QCD. The total 
CPU time consumption amounts (very roughly) 
3 • 10 8 seconds (about 10 years) if normalized to 
the use of a stand-alone 1.5 GH PC. Due to the 
heavy use of the SGI cluster (of 32 parallel SMP 
CPU of 1.5 GH frequence each) the calculation 
took about 15 calendar months. 

The corresponding calculations for the vector 
correlator will increase this demand by another 
factor of 3 (optimistically) or 10 (pessimistically). 
Clearly the combined experience of improved pro- 
grams and better hardware will lead to the de- 
sired result within the next few years. 
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